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Abstract Superhalo electrons appear to be continuously present in the interplanetary 
medium, even at very quiet times, with a power-law spectrum at energies above ~2 keV. 
Here we numerically investigate the generation of superhalo electrons by magnetic reconnec- 
tion in the solar wind source region, using the MHD and test particle simulations for both 
single X-line reconnection and multiple X-line reconnection. We find that the direct current 
electric field, produced in the magnetic reconnection region, can accelerate electrons from an 
initial thermal energy of T ~ 10° K up to hundreds of keV. After acceleration, some of the 
accelerated electrons, together with the nascent solar wind flow driven by the reconnection, 
propagate upwards along the newly-opened magnetic field lines into the interplanetary space, 
while the rest move downwards into the lower atmosphere. Similar to the observed superhalo 
electrons at 1 AU, the flux of the upward-traveling accelerated electrons versus energy dis- 
plays a power-law distribution at ~ 2 — 100 keV, f(E) ~ E~°, with a ô of ~ 1.5 — 2.4. For 
single (multiple) X-line reconnection, the spectrum becomes harder (softer) as the anomalous 
resistivity parameter a (uniform resistivity 77) increases. These modeling results suggest that 


the acceleration in the solar wind source region may contribute to superhalo electrons. 
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1 INTRODUCTION 


Electron measurements from the 3D Plasma and Energetic Particle instrument on the WIND spacecraft near 
1 AU find a superthermal component of solar wind electron population at energies above ~2 keV, denoted 


the “‘superhalo”, with a power-law (f(E) ~ E~? with 6 ~ 2.5) spectrum extending to > 100 keV, and a 
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nearly isotropic angular distribution trinli997 1998). Superhalo electrons appear to be the electron coun- 
terpart of the power-law-tail suprathermal ions above solar wind and pickup ion energies that are observed 


throughout the heliosphere at all the times (e.g.,|Gloeckler et_al.[2008;|Wimmer-Schweingruber et al.|2013). 
Using high sensitivity measurements from the SupraThermal Electron instrument (Lin et alj/2008) on the 
eet nearer boi} reported that the power-law spectral index 6 of superhalo electrons observed 


during quiet-times near solar minimum ranges from ~ 1.5 to ~ 3.4, with an average of ~ 2.35 + 0.45. The 


observed density of superhalo electrons, about 107° — 107° of the solar wind proton density, decreases with 
the decay of solar cycle, while 6 has no solar-cycle variation. Since these superhalo electrons are present 
even in the absence of any solar activity (e.g., active regions, flares, vies Wang etal suggested that 
superhalo electrons may be generated by wave-particle interactions in the interplanetary medium (IPM), or 


by nonthermal processes related to the acceleration of the solar wind. 


With a weak turbulence re | boi) proposed that superhalo electrons are accelerated 


by local resonant interactions with electron beam-excited Langmuir waves and that the dynamic equilibrium 
between these electrons and Langmuir waves predicts a power-law spectrum of E~?-3, consistent with 
observations. However, (2008) argued that the observed energy density of the Langmuir waves at 


1AU is too small to accelerate superhalo electrons by the time the solar wind reaches 1 AU. 


If superhalo electrons arise from the Sun, one possible accelerator is the magnetic reconnection in 


the solar wind source region. Many researches have examined the charged-particle acceleration by m 


Browning & Vekstei bool i | 
Turkmani et all 200€ (Carin et all 


wetall 2o10} loka etal 2011 2012 2013 


2014) |Soeised was the first to analyze the direct current (DC) electric field acceler- 
ation in the single reconnecting current sheet (RCS), by analytically solving the particle motion equations 
in the geomagnetic tail. Based on simulations with the single RCS in solar cea obtained 
a power-law energy spectrum with the index of 2.0-2.2 for the accelerated protons, consistent with the the- 


oretical prediction sy lente (1980); Wood & Neukirch| (2005) got a power-law energy spectrum with 


the index ~ 1.5 for the accelerated electrons in the reconnection outflow region. Moreover, 
boosh and|Cargill et al. (006) found a power-law distribution for the accelerated ions and electrons, from 
simulations with the stochastic development of transient RCSs. [Onofri et al. 2006) also suggested that 


the fragmented RCSs can be very efficient electron accelerator. Based on two-dimensional particle-in-cell 


(PIC) simulations of multi-island reconnection, 2010) and l2010 proposed that 


electron and ion acceleration is dominated by Fermi reflection in contracting and merging magnetic islands. 


However, the above models related to solar flares may not conform well to the superhalo electrons that are 


not correlated with flares. Due to their continuous presence in the PM, similar to the solar wind, these super- 
halo electrons could be produced by the magnetic reconnection related to the solar wind origin (Wang et al 
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If the superhalo electrons observed at 1 AU originate from the solar wind source region, then a self- 
consistent acceleration model would also include the particle-escape process from the Sun to the IPM. 
One possibility is to involve open magnetic field lines out into the IPM, as suggested by many mod- 


els of solar energetic particles from the acceleration in transient events such as solar flares and coronal 


mass ejections (Vainio et al 2003:|Arzner & Vlaho Rosdahl & Galsgaard 
Baumann & Nordlund bora. Recently, we have simulated the reconnection between 


a closed loop and an open funnel in the solar wind source region, to examine the origin of the solar wind 
2013 


test particle model (Section 2.2), and simulate the flux energy spectrum and production rate of accelerated 


). In the present study, we utilize this magnetic reconnection model (Section 2.1) and a 


electrons by the reconnection electric field in the solar wind source region (Section 3), to investigate the 


solar origin of superhalo electrons. 


2 NUMERICAL METHOD 


2.1 Numerical MHD Model 


The numerical MHD model used here has been described in details in (2013). This section 


only gives the basic features and specifies parmaters/initial conditions for this study. In order to sustain 
conversation laws and correct relationships of quantities across discontinuities in simulations, the 2.5-D 
resistive MHD equations, in the Cartesian coordinates (x, y, z) with y directed vertically, are written in a 
conservational form. By adopting reference values of the plasma density po (= 2 x 1071? kg m~3), length 


Lo (= 1 Mm), and temperature To (= 104 K), these equations are normalized as follows: 


2 Y. pu=0 (1) 
“pn +V. [puu + 10+ 5B’) — BB| = pg (2) 
Oe 1 r 
Pa ufe +p+ ZB?) - (u-B)B| = pu-g +V: (B x ni) byt Vat H+ Cy (3) 
“3 LV. (uB - Bu) = 1V°B (4) 
where 
e= sou + P+ SBP, I=V <8 (5) 


correspond to the total energy density and current density, respectively. Here, p is the mass density; u = 
(vz, Vy, Uz) is the plasma velocity; p is the thermal pressure; B denotes the magnetic field; g (= —ge,, g = 
const) is the solar gravitational acceleration; y (= 5/3) is the adiabatic index; L, represents radiative losses; 
V - q gives the anisotropic thermal conduction; Cy is the Newton cooling term; H is the parameterized 
heating term; and 77 is the magnetic resistivity. Here, we assume two resistivity models of the magnetic 
reconnection: anomalous resistivity to trigger single X-line reconnection, and uniform resistivity to trigger 
multiple X-line reconnection. 

At the Sun, once the current density J is larger than the threshold of current-driven micro-instability 


(such as ion acoustic instability), large enough local diffusion can be triggered to increase_the local 
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Büchner & Elkinaj2005)). Many theoretical analysis and simulations have reported that the anomalous resis- 
tivity strongly depends on J ae Davidson & Gladdj1975; |Petkaki & Freema dae E 
: 


e.g., in a linear or power-law form. Based on the MHD simulations by 
Sato & Hayashil (197 gai {1992} adlad (2001), here we set 7) in the anomalous resistivity model 


to be a function of the current density: 


Ko 


0, J< Je 
= (6) 
mJ = Je), I> Ie 


where 7 and qa are the resistivity parameters, and Je is the current-density threshold above which the 
anomalous resistivity is triggered. As suggested bylTreumand oot) E am boos), No and 
Je are set to be 0.001 LoVo and 0.5 BoLo ‘y7!, respectively, where Vo(= /RTo, R is gas constant) is the 
reference velocity, and Bo(= \/ZpoVo, u is magnetic permeability) is the reference magnetic strength. For 
a, we simulate for four different values, 0.5, 1.0 ne es | 2.0 load sanih. and 3.0 
1979) 

At the Sun, on the other hand, the collisional resistivity is generally too small (< 107 Q m) to be 
resolved under current computation ability, and to make dissipation difficult ee e In the 
uniform resistivity reconnection model, thus we simulate with an enhanced n, for three values of 8 x 10-4, 


2 x 1074 and 5 x 1075 (Biskamp & Welte 1080: [Samtaney et all2009 Onofri et al. 2004 2006). For a 


smaller uniform 7 (e.g., 1075), however, the upward-traveling electrons can not be accelerated to energies 
above 50 keV in our reconnection model. 

In both resistivity models, the simulation region spans —13 Mm < x < 0 Mm in the horizontal dimen- 
sion and 0 Mm < y < 15 Mm in the vertical dimension. This region is covered by a nonuniform grid in 


both dimensions, with a grid spacing of dx = dy = 25 km for 0 Mm < y < 4 Mm, ôx = dy = 50 km for 


4Mm < y < 6Mnm, and ôx = dy = 100 km for y > 6 Mm. We use the same boundary conditions defined 


For the initial conditions, we assume the plasma in hydrostatic equilibrium with the temperature of a 
hyperbolic tangent function, and derive the plasma pressure and density from the static equation and ideal 
gas equation Feaaocie oak The initial plasma velocity is set to be 0. The initial magnetic 
field is set to be an open funnel potential field given fae bood, plus a closed loop potential 


field generated by two infinite straight-line currents in the z—direction at (x = —10 Mm, y = —0.7 Mm) 


and (x = 10 Mm, y = —0.7 Mm). We also assume the presence of a relatively weak magnetic field 
component (Bo) in the z—direction, to allow the guiding center approximation for the particle orbits in the 
reconnection region, and the effective partcile acceleration E E Browning & Vekstein||2001 
Wood & Neukirchl2005 2012). 

2.2 Test Particle Approach 


Provided that the electron gyro-radius (gyro-period) is much smaller than the model scale length (char- 


acteristic time), the guiding center approximation is employed in the present simulation. The relativistic 


Numerical Simulation of Superhalo Electrons in the Solar Wind Source Region 5 


b (7) 


m (Yv))? m m (Yv 
u =up + -i [b x (b - V)b)] + cin [b x (V(KB))] + zeh x (b- V)uz)| 
z ED x (us: V)b)] + lb x (ue: Vue) (8) 
L = 1g .b— Aib -V(«B)) + (Yu )ug - ((b - V)b) + Tug: ((ug- V)b) (9) 
I ee (10) 


In the above equations, r is the electron position vector and b(= B/B) is the unit vector of magnetic 
field B; m and q are the electron mass and charge; c is the speed of light; u1 is the electron velocity compo- 
nent perpendicular to B, including E x B drifts, curvature drifts and gradient drifts; vy is the electron veloc- 
ity component parallel to B; 4p is the particle magnetic moment. T(= 1/ J1— 2/2 ) is the relativistic 
factor, where v is the electron speed, and the coefficient « equals to Vi- ur? e, where ug = E x b. 

Since the acceleration time (~ 0.001 s) of electrons is much shorter than the characteristic timescale 
(~ 20 s) of the RCS evolution in the solar wind source region, we simulate the trajectory and velocity 
of electrons with quasi-static background field snapshots from our MHD model. We spatially bi-linearly 
interpolate these MHD model results to obtain u(r), B(r) and j(r), and then calculate E(r) via E = 
—u x B + nj. The fourth-order Runge-Kutta method is employed to numerically integrate the motion 
equations (7)-(9), where the time step (ôt) is adaptive E | We also assume that electrons are 
initially distributed uniformly in the transition region, with a Maxwellian velocity distribution of T ~ 10° 
K and no bulk velocity. The rectangle in Figures 1 and 6 denotes the injection region for test electrons. In 


each case of a = 0.5, 1.0, 2.0, 3.0 and 7 = 8 x 1074,2 x 1074,5 x 1075, about 10° electron orbits are 


calculated. 


3 SIMULATION RESULTS 
3.1 Single X-line Reconnection 


For single X-line reconnection driven by anomalous resistivity, the overall evolution of the MHD simulation 
2013). This MHD simulation is not sensitive to the values 


1979; 2001), while a influences the spectral 


shape of the accelerated electrons in the test-particle simulation. In this session, we illustrate the detailed 


is similar to that in our previous work ( 


of a, consistent with previous studies (Sato & Hayashi 


simulation results for œ = 1.0, followed by brief descriptions for the other three a values. 
Figure[I]presents a MHD snapshot for a = 1.0, with the spatial distribution of the temperature T (a), 
vertical velocity Vy (b), out-of-plane current density J, (c) and total diffusive electric field Eresi (d), at 


t = 25 minutes (the early stage of reconnection). We note that carried by the horizontal flow implemented 
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Fig. 1 Spatial distributions of the plasma temperature T (a), vertical velocity V, (b), out-of- 
plane current density J, (c) and total diffusive electric field Epes; (d) at t = 25 minutes for 
the anomalous resistivity parameter a = 1.0. In panel a and b, white streamlines show the 
magnetic field lines, and black arrows indicate the plasma velocity. In panel c (d), red, green, 
blue and yellow dots denote the initial (last) positions of electrons over a time interval of 0.1 s, 
respectively, for the final energy between 2-5 keV, 6-10 keV, 11-60 keV, and 61-200 keV. This 
MHD snapshot is used for the test particle simulation, with the black rectangle denoting the 


region where the test electrons are initially injected. 


Such reconnection produces both upward and downward outflows in the reconnection region, while the 
post-reconnection pressure gradient causes a second upward flow along the newly-opened magnetic field 
lines. In the reconnection region, the large current density J forms, and the temperature there is enhanced 
due to Joule dissipation. Once J exceeds the threshold Je, the anomalous resistivity will be switched on. 


As shown in Figure Otd), therefore, the diffusive electric field Evesi(= nJ) builds up in the reconnection 
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Fig.2 (a): the trajectories in the x — y plane of four sample electrons with the final energy of 
6 keV (purple), 48 keV (blue), 107 keV (green) and 161 keV (red), superimposed in the spatial 
distribution of current density J, for a = 1.0. The colored dots indicate the electron initial 
positions. (b): their trajectories in the z — y plane. (c-d): the electron energy and Æ) versus y, 


along the trajectories of these four electrons. 


Figure[Ī{c) and (d) display, respectively, the initial and last positions of electrons over a time interval 
of 0.1 s, for the final energy between 2-5 keV (red), 6-10 keV(green), 11-60 keV (blue) and 61-200 keV 
(yellow). In the test particle simulation, initially about 10° test electrons are distributed uniformly in a rect- 
angular region containing the reconnection site (see Figure[I{b)), with a Maxwellian velocity distribution of 
T ~ 10° K and no bulk velocity. After a time interval of ~ 0.1 s, about 5 x 104 electrons (5% of the 10° test 
electrons) pass through the reconnection region, where the diffusive electric field component F (parallel to 
B) is large, and all of them are strongly accelerated to energies above 2 keV by Æ. The closer the electrons 
reach to the region with larger Ej, the higher energy they could gain (see Figure Otc) and B). After the 
acceleration, these accelerated electrons drift out of the reconnection region along the magnetic field lines 
(see Figure[Id)). About half of the accelerated electrons move upwards along the newly-opened magnetic 
field lines into IPM, together with the nascent solar wind flow driven by the reconnection (see Figure[I|d)). 
On the other hand, the other half of accelerated electrons move downwards into the lower atmosphere, and 
they would collide with the ambient dense plasma to emit Hard X-rays via non-thermal bremsstrahlung. 
Such Hard ‘bo, .., would be very weak, probably contributing to the quiet-Sun Hard X-rays (e.g., 


2010), since the observed flux of superhalo electrons at 1 AU is only ~ 10° — 10° of the peak 


flux of typical solar energetic electron events associated with Hard X-ray bursts. 


Figure[2{a) and (b) show the trajectories of four sample electrons with final energy of 6 keV (purple), 
48 keV (blue), 107 keV (green) and 161 keV (red). Although the simulation itself is two dimensional in x 
and y, the z-displacement of electrons is calculated by integrating vz over time. Figure[2{c) and (d) show, 
respectively, the electron energy and parallel electric field £ along these four trajectories. As electrons 
move along the magnetic field lines and approach the reconnection region with non-trivial Æj, they start to 


be energized. The electrons reaching the very center of the reconnection region can be accelerated (by large 
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Flux [s7 m>? keV"'] 
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Fig.3 Flux versus energy spectrum of upward-traveling electrons with the final energy from 


~ 2 keV to 200 keV for a = 1.0. The squares and asterisks represent the observed superhalo 
electrons at 1 AU from Wang et all boi), with the flux shifted by five orders of magnitude to 
compare with the simulation results. The red-dash straight lines represent a power-law fit to the 


simulations results and observations. 


the negative z-direction and reach a maximum displacement up to ~ 0.7 Mm (see Figure[2{b)). Afterwards, 
electrons leave the reconnection region and move along magnetic field lines mainly in the x — y plane; the 


acceleration dramatically decreases. 


Figures [I] and [2]suggest that the electric and magnetic field configurations built up by the magnetic re- 
connection in the solar wind source region are capable of accelerating electrons from thermal to superther- 
mal energies. Assuming a continuous electron injection into the magnetic reconnection region, we calculate 
the flux of upward-traveling electrons around x = —5 Mm and y = 10 Mm. For a = 1.0 (see Figure), the 


flux versus energy spectrum of electrons at ~ 2-100 keV fits to a power-law distribution, f(E) ~ E~?°. 


This spectral index 6 of 2.0 is consistent with the average index (2.35 + 0.45) of superhalo electrons ob- 
served in situ during quiet-time periods {wang et at 2019. 


the magnetic reconnection evolves, the spectral shape of accelerated electrons doesn’t change very much. 


We also note that in the present simulation, as 


In the anomalous resistivity model, the energy spectral shape is the most sensitive to the resistivity 
parameter a. Here, we also simulate for œ = 0.5, 2.0 and 3.0, with the other parameters fixed. Figure [4] 
compares the spatial distribution of total diffusive electric field Epes; around the reconnection region for 
a = 0.5,1.0, 2.0 and 3.0. In the reconnection with asymmetric inflow, the strong diffusive electric field 
Eyres; is built up on the strong field side of RCS. As a increases, the RCS becomes less flat and Eresi 
becomes larger, so more electrons can be accelerated to higher energies and the resultant spectrum would 
become harder (see Figure{a)). This is consistent with the simulation ie Cees ad, 


In the four a cases, the flux energy spectra of accelerated electrons above 2 keV all exhibit a single power- 
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Fig.4 A zoomed-in view of spatial distributions of total diffusive electric field Epesi at t = 25 
minutes for four simulation cases with a=0.5, 1.0, 2.0 and 3.0, in the area around the reconnection 


region. Streamlines show the magnetic field lines. 


superhalo electrons, while when a increases to 3.0, 6 decreases to ~ 1.3, harder than the observations (see 


Figure 5fb)). 


3.2 Multiple X-line Reconnection 


For multiple X-line reconnection driven by small uniform resistivity, we illustrate the detailed simulation 
results for 7 = 2 x 1074, followed by brief descriptions for 7 = 8 x 1074 and 5 x 107°. Figure[6] shows 
the modeled spatial distributions of the plasma vertical velocity V, at t = 25 minutes (the early stage 
of reconnection) for 7 = 2 x 1074. Here the RCS becomes unstable due to tearing instabilities and is 
fragmented into several magnetic islands. Such magnetic reconnection is no longer the standard Sweet- 
Parker like, and is inherently time-dependent. In Figure 6, the upward plasma outflow comes not only from 
the reconnection region, but also from the high-pressure leg of the newly-opened loops, similar to single 
X-line reconnection. Therefore, this bursty reconnection would not change the mass load to the nascent 
solar wind outflow as stated A ener mes Since the acceleration time of electrons is much smaller 
than the characteristic timescale of RCS evolution, we can still use the MHD background field snapshot to 
conduct test particle simulation, despite the time-dependent reconnection. 

In the test particle simulation, initially about 10° test electrons are distributed uniformly in a rectangular 


region denoted in Figure|6] with a Maxwellian velocity distribution of T ~ 10° K and no bulk velocity. After 
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Fig.5 a): Comparison of the flux versus energy spectra of accelerated electrons between a=0.5, 
1.0, 2.0 and 3.0, with the other parameters fixed. (b): Histogram of the observed power-law index 
of superhalo electrons Hii aa l2012 , with the colored arrows indicating the simulated 


power-law indexes for the four different a. 


are accelerated to energies above 2 keV by large Fj. After the acceleration, about half of these high-energy 


electrons move upwards along the newly-opened magnetic field lines into IPM. 


Figure[7ķa) and (b) show the trajectories of three sample accelerated electrons with the final energy of 
6 keV (blue), 30 keV (green) and 110 keV (red), respectively, in the x-y and z-y plane, for 7 = 2 x 1074. 
Figure[/[c) and (d) display the electron energy and parallel electric field Ej» respectively, along these three 
trajectories. Compared with single X-line reconnection, the trajectory and energization process of electrons 
become more complicated in multiple X-line reconnection. For the electrons trapped in magnetic islands 
(e.g., see the green and red trajectories), they gain high energy as they circle around magnetic islands and 
experience non-trivial £. Note that £ produced in multiple X-line reconnection is much ( > 10 times) 
weaker than in single X-line reconnection, due to small 7. Thus, electrons would travel a longer distance 
in the z-direction to acquire a large amount of energy. For the electrons moving freely in open field lines 
(e.g., see the blue curve), they can undergo multiple accelerations by Æ as they pass by multiple X-lines. 


However, the final energy gained by these freely moving electrons is limited due to weak Ej. 


Figure [8] compares the flux versus energy spectra of accelerated electrons between n = 8 x 1074, 
2x 1074, and 5 x 1075, with the other parameters fixed. In the three cases, the electron energy spectra above 
~ 2 keV generally fit to a single power-law function, with a spectral index ô occurring within the range of the 
observed superhalo electron indexes during quiet-time periods. As 7 decreases, the RCS becomes thinner 
and the Eresi away from X-lines gets smaller, so less electrons are produced at energies below ~ 20 keV 
(although the higher-energy electrons remain essentially unchanged) and the resultant spectrum becomes 
harder. Our simulations also show that when n decreases below ~ 10~°, the upward-traveling electrons can 


not be efficiently accelerated to energies above 50 keV, and thus the simulated energy spectrum is no longer 
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Fig.6 Spatial distributions of the plasma vertical velocity V, at t = 25 minutes for uniform 
resistivity n = 2 x 1074, with streamlines showing the magnetic field lines. This MHD snapshot 
is used for the test particle simulation, with the black rectangle denoting the region where test 


electrons are initially injected. 


4 SUMMARY AND DISCUSSION 


In this study, we investigate the generation of superhalo electrons under the magnetic reconnection model 
for the solar wind origin. Using self-consistent electric and magnetic fields obtained from the MHD re- 
connection simulation, we conduct test particle simulations to study the acceleration of electrons in solar 
wind source region, for both single X-line reconnection driven by anomalous resistivity and multiple X-line 
reconnection driven by small uniform resistivity. We find that the superhalo electrons may be contributed 
by the DC electric field in the magnetic reconnection in the solar wind source region. 

The simulation results show that in both reconnection models, electrons with an initial Maxwellian 
velocity distribution of ~ 10° K can be accelerated to high energies, up to hundreds of keV, mainly by 
the DC electric field in the magnetic reconnection. For single X-line reconnection, electrons gain higher 
energy as they get closer to the reconnection region, while for multiple X-line reconnection, electrons gain 
high energy as they are trapped in and circle around magnetic islands. Afterwards, the accelerated electrons 
follow magnetic field lines to drift out of the reconnection region. About half of the accelerated electrons 
propagate outwards along the newly-open magnetic field lines into IPM, together with the nascent solar 
wind flow driven by the reconnection, while the other half move downwards into the lower atmosphere. 

In both reconnection models, the energy spectrum of the upward-traveling electrons generally fits well 


to a single power-law at energies of ~ 2 — 100 keV. For single X-line reconnection, the simulated spectral 
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Fig.7 (a): the trajectories in the x — y plane of three sample accelerated electrons with the 
final energy of 6 keV (blue), 30 keV (green) and 110 keV (red), superimposed in the spatial 
distribution of current density J,, for uniform resistivity 7 = 2 x 1074. The colored dots indicate 


the electron initial positions. (b): their trajectories in the z — y plane. (c-d): the electron energy 


and F versus y, along the trajectories of these three electrons. 
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Fig.8 Same as Figure[5]but for uniform resistivity 7=8 x 1074, 2 x 1074, and 5 x 1075. 


served during quiet-time periods (wang ct atlloor2). For multiple X-line reconnection, the simulated 6 is 


~ 1.5 — 2.0, within the index range of superhalo electron observations. Among the model parameters, the 
resistivity can significantly affect the RCS that, in turn, changes the energy spectrum of accelerated elec- 
trons. For single X-line reconnection, as the anomalous resistivity parameter œ increases, more electrons 
are accelerated to high energies, and thus the simulated spectrum becomes harder. For multiple X-line re- 
connection, as uniform resistivity 7 increases, more electrons can be accelerated to energies below 20 keV 
and the simulated spectrum becomes softer. 

Since test particle simulation is not self-consistent, we cannot directly estimate the number density of 


accelerated electrons and thus their density ratio to the solar wind density in the solar wind source region. 
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only ~10% of the electrons passing through the magnetic reconnection region can be accelerated to energies 
above 2 keV. Using the physical parameters of the MHD simulation as well as this assumption, the density 
of upward-traveling accelerated electrons above 2 keV is estimated to be ~10~ 1° of the nascent solar wind 
flow at y= 10 Mm (see Appendix for details). This is close to the in situ superhalo electron observations 
len if the number ratio of superhalo electrons to solar wind plasma doesnt vary significantly 
en route to 1 AU. But note that it is unknown whether the 10% acceleration efficiency assumption is valid 


in the solar wind source region. 


We should note that in the present MHD simulation proposed by [Yang et all (2013), only one closed 


loop is considered to reconnect with open field lines in the funnel rooted at a chromospheric network con- 
junction. After the magnetic flux in this closed loop is used up, the reconnection process will cease. To 


form a continuous and relatively uniform solar wind as observed in the IPM, we need to consider a large 


number of independent micro reconnection events that happen in a limited time period and in many funnels 
(Tu et a. 


2005tb, 


this scenario assumes that the successive impulsive reconnection events take place at the interfaces between 


in the wind source region. Like the reconnection scenario in the solar wind source region 


meso-scale closed loops within the intra-networks and open funnels rooted from the network junctions, to 


account for the continuous presence of superhalo electro population in the IPM. Those successive impul- 


sive reconnection events may cause local transient events such as bi-directional plasma jets and network 
brightening 2000), as well as large-scale quasi-steady Doppler blue shift in 
the higher transition region and corona 2005 at 2010). In the future, we 


will develop a 3-D simulation model to describe the formation of continuous wind from many intermittent 


micro reconnections and accordingly develop a more realistic particle acceleration model. 


This study is mainly aimed at investigating whether the superhalo electrons in the IPM could originate 
from the Sun. The present model, although simple, is the first to numerically validate this scenario. In 
future, we can improve this model by considering a comprehensive acceleration, instead of only the DC 
electric field acceleration. Other mechanisms, such as the turbulence (or stochastic) acceleration, shock 
acceleration, or collapsing magnetic trap acceleration, could also take place in the reconnection region 
eect, Including these mechanisms would likely allow more electrons to be accelerated 
to higher energies, increasing the density ratio between the accelerated high-energy electrons and solar 
wind. Moreover, the Fermi acceleration by reflection from contracting and merging magnetic islands could 
occur in the multiple X-line reconnection trace O N As a result of the mismatch between the 
spatial and temporal scales of the MHD fields and those of electron motion, our MHD models are unable 


to address this Fermi acceleration; instead, the PIC simulation is usually used to study it beaa 2010 
2010). At the solar wind source region, however, the very small ion inertial length makes the PIC 


simulation subject to resolution constraints, reducing the possible physical box size that can be simulated 
to far below the length scale of the reconnection region. In future, we may also include PIC simulation to 
examine the influences of the Fermi acceleration by contracting and merging magnetic islands on the results 


presented here. 


Also, we can improve the model by taking into account the effects of Coulomb collisions at the Sun 
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model would require simulations to be carried out on longer time scales. In simulations with collisions, 
have found that the effect of collisions becomes dominant with time, since the 
source of acceleration (strong electric field) is transient and thus gradually disappears, while the source 
of energy losses (coulomb collisions) is always present. Therefore, at the early stage of reconnection, the 
energy spectra are rather similar to those obtained in simulations with no collisions, while at the later stage 
when electric fields are gradually decaying, the collisions become dominant. In our test particle simulation 
for the early stage of reconnection occurring at transition region, we will check the effects of Coulomb 
collisions on the acceleration of electrons, although such effects may be insignificant. In addition, during 
the interplanetary propagation, reflection by the interplanetary shocks (e.g., CIR shocks), and/or scattering 
by wave-particle interaction (e. o konalo Mocks e alpoos can isotropize the angular distribu- 
tion of superhalo electrons, to form a nearly isotropic distribution observed at 1 AU. Also note that the 
evaluated flux and density of the accelerated electrons in the present model are based upon the value of ac- 
celeration efficiency, which can be sensitive to various MHD model parameters, especially to the resistivity 
parameters. In future, we will consider these respects, to compare the simulations and the observations in 


details. 
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Appendix A: APPENDIX 
A.1. Single X-line reconnection 


In the test particle simulation, initially about 10° test electrons are distributed uniformly in a rectangular 
region containing the reconnection site (see Figure[I|b)), with a Maxwellian velocity distribution of T ~ 10° 
K. After a time interval of ~ 0.1 s, about 5 x 104 electrons (5% of the 10° test electrons) pass through the 
reconnection region, where the diffusive electric field component F) (parallel to B) is large, and all of them 
are strongly accelerated to energies above 2 keV by Ej. After the acceleration, these accelerated electrons 
drift out of the reconnection region along the magnetic field lines (see Figure[I{d)), with an average final 
energy Kave1 ~ 4 keV and average final velocity Vave1 ~ 3.5 x 107 ms~!. About half of the 5 x 104 
electrons move upwards along the newly-opened magnetic field lines into IPM, together with the nascent 
solar wind flow driven by the reconnection (see Figure[I{d)). 

Using the physical parameters of the MHD simulation, we can estimate the actual flux, density and total 
energy gain of the accelerated electrons above 2 keV. In the simulated rectangular region (see Figure 1(b)) 
with an area Sy = 4.5 Mm x 3.5 Mm in the x — y plane and a depth lı = 1 cm in the z—direction, the 


actual total number of initial thermal electrons is 
Ny = no x Sı x l = 1.6 x 107°, (A.1) 


where no = 10° cm~® is the background plasma density. According to the test particle simulation, about 
5% of the total electrons would pass through the reconnection region and be all accelerated within a time 


interval of ~ 0.1 s. In a self-consistent simulation, however, not all of them would be accelerated to high 
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~ 10% of the electrons passing through the magnetic reconnection region can be accelerated to energies 
above 2 keV. In the present simulation, therefore, the physical production rate of the > 2 keV electrons can 
be estimated as: 


Pa = (10% x 5% x N1)/0.1 s = 8 x 10” 571. (A.2) 


According to the test-particle simulation, about half of the accelerated > 2 keV electrons would move 
upwards along the newly-opened magnetic field lines into IPM. At y = 10 Mm, these electrons cross an 
area Ssupı =~ 0.2 Mm (defined as the width at the 1/10-fold peak intensity of electron spatial distribution 
at the x—direction) x 1 Mm (the 1/10-folding width in the z—direction). Thus, the flux of upward-traveling 


> 2 keV electrons is fsupı = 0.5 x Pei /Ssup1, and their number density is 


Nsupl = fsup1/Vavei =~ 5.7 X 104 m`’, (A.3) 


At y = 10 Mm, the average number density of simulated solar wind outflow is n5w1 =~ 2.7 x 1014 m~’ 


from the MHD simulation. Then the simulated ratio nsupı /Nswiis 2x 10~!° at this altitude. 


A.2. Multiple X-line reconnection 


In the test particle simulation, initially about 10° test electrons are distributed uniformly in a rectangular 
region denoted in Figure [6] with a Maxwellian distribution of T ~ 10° K. After a time interval of ~ 0.1 
s, about 3 x 10* electrons (3% of the 10° test electrons) pass through the reconnection region, where the 
diffusive electric field component Æj (parallel to B) is large, and all of them are strongly accelerated to 
energies above 2 keV by F. After the acceleration, these accelerated electrons drift out of the reconnection 
region along the magnetic field lines, with an average final energy Kave2 ~ 3 keV and average final velocity 
Vave2 ~ 3.2 x 107 m s7}. About half of the 3 x 104 electrons move upwards along the newly-opened 
magnetic field lines into IPM, together with the nascent solar wind flow driven by the reconnection. 

Using the physical parameters of the MHD simulation, we can estimate the actual flux, density and total 
energy gain of the accelerated electrons above 2 keV. In the simulated rectangular region (see Figure|6) with 
an area S2 = 3 Mm x 2 Mm in the x — y plane and a depth l2 = 1 cm in the z-direction, the actual total 


number of initial thermal electrons is 
No = no Xx Sq X lg = 6 x 1074, (A.4) 


where no = 108 cm~® is the background plasma density. According to the test particle simulation, about 
3% of the total electrons would pass through the reconnection region and be all accelerated within a time 
interval of ~ 0.1 s. In a self-consistent simulation, however, not all of them would be accelerated to high 


energies. Then the physical production rate of the > 2 keV electrons can be estimated as: 
Pez = (10% x 3% x N2)/0.1 s = 1.8 x 10% s7}, (A.5) 


for a 10 % acceleration efficiency in the number of electrons passing through the reconnection region 


(Baumann & Nordlund2012). 


According to the test-particle simulation, about half of the accelerated > 2 keV electrons would move 


upwards along the newly-opened magnetic field lines into IPM. At y = 10 Mm, these electrons cross an 
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at the x—direction, from the test particle simulation) x 0.2 Mm (the 1/10-folding width in the z—direction). 
Thus, the flux of upward-traveling > 2 keV electrons is fsup2 = 0.5 x Peo if Ssup2, and their number density 
is 

Nsup2 = fsup2/Vave2 =~ 3.4 x 10* m™, (A-6) 


At y = 10 Mm, the average number density of simulated solar wind outflow is nsw2 =~ 3.4 x 1014 m~’ 


from the MHD simulation. Then the simulated ratio nsup2 /Nsw2 IS 10-19 at this altitude. 


